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ABSTRACT 

Cusp-core problem is a controversial problem on galactic dark matter haloes. Cosmo- 
logical A^-body simulations has demonstrated that galactic dark matter haloes have 
a cuspy density profile at the centre. However, baryonic physics may affect the dark 
matter density profile. For example, it was suggested that adiabatic contraction of 
baryonic gas makes the dark matter cusp steeper. However, it is still an open question 
if the gas falls into the galactic centre in smooth adiabatic manner. Recent numeri- 
cal studies suggested that disc galaxies might experience clumpy phase in their early 
stage of the disc formation, which could also explain clump clusters and chain galaxies 
observed in high redshift Universe. In this letter, using numerical simulations with an 
isolated model, we study how the dark matter halo responds to these clumpy nature 
of baryon component in disc galaxy formation through the clump cluster phase. Our 
simulation demonstrates that such clumpy phase leads to a shallower density profile 
of the dark matter halo in the central region while clumps fall into the centre due 
to dynamical friction. This mechanism helps to make the central dark matter density 
profile shallower in the galaxies whose virial mass is as large as 5.0 x 10^^ M©. This 
phenomenon is caused by reaction to dynamical friction of the stellar clumps against 
the dark matter halo. The halo draws the clumps into the galactic centre, while kine- 
matically heated by the clumps. We additionally run a dark matter only simulation 
excluding baryonic component and confirm that the resultant shallower density profile 
is not due to numerical artifact in the simulation, such as two-body relaxation. 

Key words: methods: numerical - galaxies: formation - galaxies: evolution - galaxies; 
haloes - galaxies: spiral. 



1 INTRODUCTION 



Cosmological Ai'-body (dark matter only) simulations with 
a large number of particles have been carried and demon- 
strated that the density profiles of the dark matter (DM) 
halo have cusps at their central region (e.g. Navarro et al.l 
19971: iFukushige fc Makinolll99l l200ll: iKlvpin et al.ll200ll: 



Navarro et all 12004 ISpringel et all l2008l : iNavarro et all 



2010l : llshivama et al.l l201lh . However, it has been obser- 
vationally suggested that DM haloes of dwarf and low 
surface brightness (LSB) galaxies seem to have a central 
constant density region (cored ha lo), not a central cusp 
("Gilmorc ct al. H2007l : Ide Blok|[2010l . and references therein; 
foh c t al. 20lJ). 

In the real Universe, the galactic DM haloes might 



be affected by baryonic physics. Recently, some studies 
suggested that a probable mechanism to turn the cusp 



bursty star formation in dwarf galaxies ( 


Navarro et al.lll996l: 


Mo & Mad 20041: Read & Gilmord 2005 


; Mashchenko et al. 


20081: 


Governato et all 201C 


: iPasetto et all l2010l: lOh et al. 


201 J; 


Pontzen & Governato 


I2OII). Others oroDOsed core 



creation by angular momentum tran sfer from gas c louds 



globular clusters , dwarf galaxies 



iTonini et aLll2006l:lGoerdt et al 



201c 



El-Zant et al] 



Cole et all 



2011 



200 ll : 



) and 



galac tic bar (|Weinberg fc Kat j |2002| : lEl-Zant fc Shlosmanl 
|2002| ). These baryonic structures lose angular momentum 
by dynamical friction against DM halo. Then, these studies 
suggested that the central DM cusp is dynamically heated 
and smoothed out by reaction of the dynamical friction. 
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However, most of these discussions of core creation were 
concentrated on less massive galaxies of which masses are 
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< 10^~^° M0. It is because, unlike massive systems, rela- 
tively shallow potential wells of these systems are suscep- 
tive to dynamica l perturbations (|Mac Low fc Ferraralll999l : 
iDel Popoldboill ). 

Since the systems of dwarf and LSB galaxies have a high 
mass-to-light ratio even in the central region, the DM pro- 
file is more promisingly measurable by kinematics of baryon 
tracers. On the other hand, most galaxies are dominated 
by baryon in their central region. Therefore, it is difhcult 
to clarify DM density profile by any ob servation. Ifowever , 
although there have been fewer studies, ISpano et al.l (|2008l ) 
has discussed from observation of rotation velocity curves 
that high surface brightness disc galaxies whose types range 
from Sab to Im also favors the cored DM halo rather than 
the cuspy density profile, as well as dw a rf and LSB galaxies. 
Contrary to the result of ISpano et al.l (|2008l 'l. on the other 
hand, some theoretical studies have proposed that central 
density profile of DM halo of massive galaxies become steeper 
rather t han shallower by adiabatic contraction of bary onic 
gas (e.g. lBlumenthal et al.lll986l : ICnedin fc Zhadl2002l ). 

The theory of adiabatic contraction assumes slow and 
smooth gas infall. However, some galaxies observed at 
high redshift Universe have displayed clumpy stellar struc- 
tures in their system, which ar e referred to as clump clus- 
ters and chain galaxies (e. g. Ivan den Bergh et alT 1 19961 : 



lElmegreen et al.l |2004 l2009l : lGenzel et al.ll201ll '). Numeri- 
cal simulations suggested that since the highly gas-rich disc 
is expected in early stage of the disc galaxy formation, 
the clumpy str uctures form due to the instability of a c- 
creting gas (e.g. lNoguchilll998l.[l999l:llmmeli et al.H2004allbl: 



Bourna ud et all l2007l: lAgertz et~aLl |2009| : ICeverino et all 
'2OIQ . i201ll : lAumer et al.lbOld '). Thev suggested that these 



clumpy galaxies would evolved into current disc galaxies. 
In this scenario, the stellar clump has muc h larger mass (~ 
10^~ ^ M.q) than ordinary star clusters fe.g. lElmegreen et al.l 
I2OO9I ). Therefore, the clumps are expected to fall into the 
galactic centre by dynamical friction. This clumpy manner 
of baryon accretion may be far from adiabatic. Therefore, it 
is an open question how the DM responses to such clumpy 
disc formation. 

Although some studies mentioned above have investi- 
gated the core creation due to angular momentum transfer 
in various cases of dynamical heating sources, it is yet to be 
studied how the DM halo reacts to the disc galaxy formation 
through clump cluster phase. In this study, we run a numer- 
ical simulation for disc galaxy formation in an isolated DM 
halo and study how DM halo cusp would be affected in disc 
formation through the clump cluster phase. 



2 SIMULATION 

We use an Ai'-bo dy/Smoothed Par t icle H ydrodynamics 
(SPH) code, ASURA (Saitoh et al.ll2008l . [2OO9I '). Gravitational 
interaction is calculated b y parallel Tre e with GRAPE 
(GRAvity PipE) method (jMakind |2004 ). where we use 
Phantom-GRAPE (Nitadori et al . in prep.). Gas dynamics 
is ha ndled by the standard SPH (|Monaghanlll992l : ISpringell 
I2OIOI ). The number of neighbouring particles of each SPH 
particle is set to 32±2. Equations of motion of particles 
are integrated by the leap-frog method with the individ- 
ual time-step scheme. We here adopt the time-step lim- 



iter (|Saitoh fc Makind |2009| ) which enforces the neighbor- 
ing particles to have similar time-steps so that SPH with 
individual time-steps can solve the evolution of the gas 
in the s hocked regions correct ly. In addition, the FAST 
method (|Saitoh fc Makindl20ld ') is also adopted in the time- 
integration; this algorithm accelerates the simulations of 
self-gravitating gas by integrating Hamiltonian split into the 
gravitational potential and others (kinetic and internal en- 
ergies) with different time-steps. 

Radiative cooling of gas is solved by assuming an op- 
tically thin cooling function depending on metallicity of 
gas particles, which covers a wide temperature range of 
10 - 10** K (|Wada et al.ll2009l ). Radiative heating due to 
far-ultraviolet radiation and type-II SN feedback are also 
taken into account. Un l ike fe e dback impleme n tation s of e.g. 
iThacker fc CouchmanI (|200ll ): IStinson et al] (|2006l ). we do 
not halt the radiative cooling in thermalized regions due 
to SNe. In modeling star formation and the SN feedback, 
we adopt a simple stellar population approximation with 
the Salpeter initial mass function (IMF) (|Salpctcr,,1955i ) in 
which stellar mass ranges from 0.1 M© to 100 M©. A gas par- 
ticle spawns a stellar particle whose mass is set to have 1/3 
of the original gas particle under the criteria: 1) pgas >100 
aim/cm^, 2) temperature T < 100 K, 3) V-v< and 4) 
there is no SN around the particle. The loca l star formatio n 
rate is assumed to follow the Schmidt law l|Schmidtlll959l ). 
with the dimensionless star formation efficiency C* = 0.033. 
Since in our simulation we assume that stars form from quite 
dense and cold gas, the global properties of star formation 
such as star format i on rate or the S chmidt-Kennicutt rela- 
tion (|Schmidtlll959l : iKennicut^l 19891) i s fairly insen sitive to 
the adopted value of C, (jSaitoh et aLllioOsi . |2009| ). We as- 
sume that stars heavier than 8 M© cause type-II SNe and 
thermally heat up ambient gas. 

Our initial condition follows the spherical model of 
iKaufmann et aL 1 (Hope, ,2007,) that was used to study the 
formation of disc galaxies in an isolated environment. We as- 
sum e an equilibrium sys tem with Navarro- Frenk- White pro- 
file (|Navarro et al.|[l997l . hereafter NEW) with a virial mass 
Mvir = 5.0 X 10" Mq, a virial radius Tvir = 1.67 x 10^ kpc 
and a concentration parameter, c = 6.0. Baryon mass frac- 
tion of the system is set to 6 per cent. At the initial phase, 
all of baryon is in the gas phase. Density profiles of DM 
and gas components are initially the same. The gas is ini- 
tially primordial (zero-metal gas) and assumed to have virial 
temperature. The gas sphere has a rotation fo llowing a spe- 
cific a ngular momentum distribution of j cx r ([Bullock et al.l 
I2OOII ) normalized by a spin parameter of A = 0.1, while the 
DM halo does not rotate. The DM halo is represented by 10^ 
coUisionless particles with a softening length of eoM = 8 pc, 
whereas the gas in the halo is expressed by 5.0 x 10® SPH 
particles with egas ~ 2 pc. Star particles formed from gas 
particles also have the same softening length as the gas par- 
ticle. 

In addition to the fiducial run, we perform a pure TV- 
body run in order to make sure that change in the central 
density profile is not due to numerical artifact, such as two- 
body relaxation. The A'^-body run uses the same density 
profile as that of the fiducial run, which is represented by 
only the DM particles with the same resolution with no net 
rotation. 
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Figure 1. Stellar surface density maps from face-on view in central 60 X 60 kpc region. Time in unit of Gyr is indicated on right 
bottom in each panel. 



3 RESULTS 

We run the simulation for 6 Gyr and show snapshots during 
the simulation in Fig[l] The galaxy simulated settles into a 
stable state after 5 Gyr, in which the density map unchanges. 

At t < 2 Gyr, many clumps form in the disc while the 
surrounding gas is falling onto the disc. In this phase, since 
the disc is highly gas-rich, radiative cooling enhanced in the 
dense gas induces self-gravity instability and results in the 
clump formation. This clumpy phases in the initial gas rich 
disc are seen in previous numerical studies and suggested 
to explain th e clump cluste r s observed at high redsh ift 
Universe ( e.g. iNoguchil Il998l. Il999l: llmmeh et al.l l2004allbl: 



10' 



Boumaud et all l2007l: lAgertz eraP l2009l : ICeverino et all 
2010l . I2OIII : I Aumer et al.l I2OI0I ) . In our simulation model, 
the gas infall ceases around f ~ 2 Gyr. The clumps decrease 
in number because of merger with on e another and tidal 
disruption (|Shlosman fc Noguchi|[l99^ ). Finally, all clumps 
me rge into a singl e "clu mp-origin" bulge as first suggested 
by iNoguchil (| 19981 . 1 19991 '). The galactic disc in our simula- 
tion is formed from clumps disrupted by galactic tide, while 
clumps migrate to the galactic center. At t = 6 Gyr, eighty 
per cent in gaseous mass had turned into stars. There is 
little gas in halo region. 

The bulge has a mass of 4.5 x 10^ M©, a half-mass ra- 
dius of Te = 450 pc and a mean surface brightness within 
of 7e — 3.5 X 10"^ Lq/pc^ under an assumption that stellar 
mass-to- light ratio is 1.0. These values are consistent with 
an re-Je relation observed am ong bulges and elliptical galax- 
ies (e.g. iBurstein et al.lll997l ). We find the transition radius 
between the disc an d bulge components at i? ~ 1 kpc on 
the surface density. iBournaud et al.l (|2009l ) has suggested 
that the disc formed from the clumps is a thick disc. In our 
simulation, a thin disc hardly forms and a thick disc domi- 
nates over the thin disc. The thick disc has a scale height of 
1.41 kpc, a scale radius of 9.85 kpc and a nearly exponential 
density profile. We will report all the details of the bulge 
and disc formation in our simulation in forthcoming paper 
(Inoue & Saitoh, in prep). 
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Figure 2. Time-evolution of density profiles of the DM halo. 
The line labelled as 'N-body 6Gyr' indicates result of the DM 
only simulation. The centre of the halo is defined to be the centre 
of mass of the DM particles. 



In the disc formation scenario through clump cluster, 
large clumps are sucked into galactic centre by dynamical 
friction against baryon and DM particles. Therefore, the 
DM halo must be kineniatically heated by the clumps as 
reaction to the dynamical friction. We show time-evolution 
of DM density profile in this simulation in Fig[2l As seen 
in the figure, DM density profile becomes nearly constant at 
the central region, r < 3 kpc, i.e. cored. The cored region ex- 
pands until t = 4 Gyr. However, the central cusp revives at 
t = 6 Gyr due to the bulge formation: settlement of a heavy 
bulge makes potential well deeper and drags matter into the 
centre. The resultant cusp slope of DM density profile is sig- 
nificantly shallower than the initi al state and central den sity 
concentration is greatly reduced (|Elmegreen et al1l2008l ). 

We perform a pure A'^-body run additionally in order to 
discuss artificial effect such as two-body relaxation. This N- 
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body system must be in the equilibrium state as long as the 
relaxation process does not work. The result is superposed 
on Fig[2] The density profile of the A'^-body run is consis- 
tent with the initial state and unchanged even after 6 Gyr, 
demonstrating that the shallowed DM density profile in the 
fiducial run is not caused by the artificial relaxation. 



4 DISCUSSION AND SUMMARY 

In this letter, we first discussed the cusp-core problem in 
the formation stage of a disc galaxy through clump cluster. 
As seen in Fig[2] the central DM density slope is clearly 
shallowed by the stellar clumps falling into the galaxy centre. 
This would be caused by reaction to dynamical friction on 
the clumps. This mechanism is basically the same as in the 
cases o f other dynamical sourc e s discussed in the previous 
works (IWeinberg fc Kat3 12002|: lEl-Zant fc Shlosmanl [2002 : 



El-Zant et all I2OOII : iTonini et al.ll2006l : iGoerdt et al.|[201C : 
Cole et al.ll201ll ). The result suggests that DM haloes would 



be cored in the clump clusters and chain galaxies in the 
high redshift Universe. In our simulation, the DM cusp has 
been smoothed out in the early stage of their formation (the 
clump cluster phase) and a shallowed DM cusp remains at 
the later epoch. This process can be effective also in current 
disc galaxies like the Milky Way, if it used to be a clump 
cluster and if there is no extra dynamical perturbation. 

Our simulation model is missing some physical pro- 
cesses possible in cosmological context, such as galaxy 
merger and accretion flow. However, at first, recent cosmo- 
logical simulations has demonstrated that disc galaxies in 
high-redshift Universe are mainly fed by smooth a ccretion 
flow, not by galactic merger (e.g. iDekel et al.ll2009l l. Then, 
we expect that the lack of the merger events would not spoil 
reasonableness of our simulation. Next, these cosmological 
simulations showed that the smooth accretion is narrow cold 
gas stream, which is different from our simulation model: 
isolated collapse model. Nevertheless, simulations using the 
isolated collapse model can also succes sfully reproduce the 
clumpy nature of hig h -redshift galaxies ([Noguchil 19981 . 1 19991 : 
llmmeh et a"l]|2004al lbl: [Aumer et al.ll2010l '). Therefore, we ex- 
pect that the isolated collapse model would not be a seri- 
ously incorrect assumption for qualitative discussion. 

Our simulation result indicates a possibility that pro- 
genitors of current disc galaxies have a cored DM halo even 
before the disc forms. However, our simulati on model is lack- 
ing t he later accretion to form a thin disc jBournaud et al.l 
I2OO9I ). Such a late accretion may also affect the DM den- 
sity. For example, the dynamically cold thin disc may forms 
spiral ar ms and/or galactic bar, which c an form a pseudob- 
ulge fe.g. iKormendv fc Kennicuttll2003 ). It is also possible 
that the pseudobulge steepens the DM density slope. Mean- 
while, the galactic bar may also flatten the densit y slope 
l|Weinberg fc Katj|2002l : lEl-Zant fc Shlosmanl 120021 ). Then, 
we have to note that it is still a open question whether DM 
haloes of current spirals are cored or cuspy. 

The cored DM density profile brings important in- 
fluences on dynamics; rotatio n cur ve ^ orbit o f star 
clusters (iTremaine &: Weinbere 1984 : Weinberg 19891: 



Hernandez fc Gilmoreljl998l : iGoerdt et al.ll2006l : iRead et all 

effect 



20061: llnouel l2009al lbh. tidal effect on the clusters 
I Huang et aLlbOOSl ). moreover, intensity of 7-ray emission (if 



annihilating radiatio n of DM was possib le and observable). 
As mentioned in ^1. ISpano et al.l (|2008l ) has discussed that 
high surface brightness disc galaxies favour the cored DM 
halo rather than the NFW cuspy profile. The 7-ray emission 
from the Milky Way centre has no t been clearly observed yet 
by ea rly result of Fe rmi telescope (|Johannll2009l:IVitale et al] 
I2OO9I , although see iGoodenough fc Hoopeij |2009| ) although 
this only means that positive evidence of enough high den- 
sity of DM has not been observed in the Galactic Centre, 
since there are some other physics to explain the no detec- 
tion. 

In our simulation, the bulge formation reconstructs the 
DM cusp at the final state although the slope is shallower 
than the NFW cusp. This result implies that DM density 
profiles in current disc galaxies would depend on size (mass) 
of the bulge: even if clump clusters make the central DM 
density cored, formation of a bigger bulge leads to the re- 
vival of steeper central cusp. Then, if we naively assume that 
all bulges are formed by clump mergers in a clump cluster, 
early-type spirals may have a cuspy DM profile since they 
tend to have a large bulge. On the other hand, since late- 
type spirals, like LSB galaxies, tend to have a smaller bulge 
(or be bulgeless), they may have a shallower cusp (or cored) 
DM halo. However, clump clusters may inevitably lead to 
formation of the clump-origin bulge. Therefore, the galaxy 
formation scenario through clump cluster phase may not be 
eligible for f ormation of the late -type and bulgeless galax- 
ies. However. ICenzel et al.l (|201ll ) has recently observed very 
strong outfl o w fro m the stellar clumps in clump clusters. 
ICenel et al.l (|2010l ) adopted such strong outflow model in 
their simulation and suggested that the clumps would be de- 
stroyed due to strong outflow from the clumps before reach- 
ing the galactic centre. In such a case, some disc galaxies 
forming through clump cluster phase may be able to evade 
the bulge formation and the central cores may remain. 

However, the clumps in our simulation seem long-lived, 
in contrast with t h e short-lived clump s suggested by e.g. 
iGenel et al.l (|2010l ): iGenzel et all (|201ll) . The SN feedback 
model adopted in our simulation can not reproduce such de- 
structive outflow. Although to investigate what is the cause 
of the difference is beyond the scope of this letter, we spec- 
ulate that IMF may be largely different from the Salpeter 
IMF in gas-rich circumstance of the clumpy galaxies. If the 
clumpy galaxies have a top-heavy IMF, the observational re- 
sult of such powerful outflow may be explainable. Addtion- 
ally, it is also discussed that SN feedback in simulations can 
not reach very high temperat ure under the simple st ellar 
population approximation (e.g. ISaitoh fc Makindfioiol ). 
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